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We showed that in Lamb-Dicke regime and under rotating wave approxi- 
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\^ ' mation, the dynamical behavior of two trapped ions interacting with a laser 

o ■ 

beam resonant to the first red side-band of center-of-mass mode can be de- 
' scribed by Jaynes-Cummings Model. An exact analytic solution for this kind 

of Jaynes-Cummings model is presented. The results showed that quantum 
^ ' collapses and revivals for the occupation of two atoms, and squeezing for 



vibratic motion of center-of-mass mode existed in both two different types 



^ ' of initial conditions. The maximum momentum squeezing for center-of-mass 



mode in these two types of conditions are found to be 42.4% and 43.8% re- 
spectively. The coherence, in the first type of initial conditions can keeps long 
times, and in the second type of initial conditions, a concrete form of coherent 

state is obtained, when the initial average number is very small. 
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I. INTRODUCTION 



Jaynes-Cummings model(JCM) [1] originating from quantum optics that describes the 
ideahzed situation of resonant interaction of an undamped two-level system with a single, 
quantized harmonic oscillator turned out to be one of the fundamental models in NMR, 
quantum optics, quantum electronics and resonance physics. This model can be exactly 
solved and predicts some interesting nonclassical effects both with respect to the atomic 
states and the boson mode, such as atomic revivals[2] and squeezing[3]. So this model has 
been studied extensively in both theoretical and experimental field. 

In the field of cavity QED, experimental realization of JCM with dissipation has been 
demonstrated both in microwave and in optical domain. Observation of quantum revivals 
and nonclassical photon statistics with Rydberg atoms and microwave cavities has been 
reported [4], and in the optical regime vacuum Rabi splitting has been observed [5]. These 
experimental progress stimulated a manifold of further theoretical studies, for example of 
the quantum correlation between two interacting subsystems [6]. Moreover, related models 
have been rendered successively, such as nonlinearly coupled [7] and multiquantum JCMs 
[8], Raman- type [9], and three-level models [10]. 

Another fundamental system to exhibit, in appropriate limits, JCM dynamics is a 
trapped ion interacting with laser beam. The quantized center-of-mass motion of the ion in 
the trap potential plays the role of the boson mode, which is coupled via the laser to the 
internal degrees of freedom. An analogy between an undamped trapped ion and the JCM 
has been noted by Blockley, Wall, and Risken [11], and by J. I. Cirac, et al., for damped 
trapped ion [12]. A study for nonlinear multiquantum JCM dynamics of a trapped ion is 
presented by W. Vogel and R. L. de Matos Filho [13]. 

In recent experiments of trapped ions, entangled states up to four ions have been re- 
alized, and some applications, such as the proof of Bell's inequality and decoherence-free 
quantum memory, have been confirmed[14]. These improvements offered some advantageous 
conditions for studying properties of multi-ion JCM. In this paper we will demonstrate that 
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a system comprising two trapped ions interacting with laser beam resonant to the first red 
side-band of center-of-mass mode, in Lamb-Dicke regime, shows a JCM-hke dynamics. Col- 
lapses and revivals of the occupation of internal states, and squeezing of the vibratic motion 
of center-of-mass mode, can be observed in two different types of initial condition. The max- 
imum momentum squeezing in these two different types of initial condition are 42.4% and 
43.8% respectively. Coherence of center-of-mass mode is also discussed in these two types 
of condition. In the first type of condition, e.g. the internal state of the two atoms being 
in ground and the vibratic motion of center-of-mass mode being a coherent state, coherence 
can keeps long times when the initial average number is very small. In the second type of 
condition, e.g. the internal state being a superposition and motional state being a vacuum, 
a concrete form of coherent state for center-of-mass mode is obtained in another case of very 
small phonons. These properties for two trapped ions are useful for quantum computation. 
Because quantum manipulation involves at lest two-bit operations, and only we know the 
dynamical properties of the ions very clearly, can we manipulate them well and truly. In 
addition, the model for this system may be realized in present experimental conditions [14], 
which will result in directly the tests of the predictions from this paper. 

This paper is organized as follows. In section II, we firstly deduce the JCM model for 
n trapped ions interacting with lasers in appropriate conditions, and then present its exact 
solution for the case of two trapped ions. In section III, we reveal the collapses and revivals 
for the occupation of internal states of the two trapped ions in two different types of initial 
condition. The coherence and squeezing for the vibratic motion of center-of-mass mode are 
discussed in section IV. Finally a summary and some conclusions are given in section V. 

II. INTERACTION MODEL 

For the sake of generality, we consider for a while n ions trapped in a linear trap which are 
strongly bounded in the y and z directions but weakly bounded in a harmonic potential in x 
direction. Assuming that the n ions are illuminated simultaneously with a dispersive beam 



3 



(or beams) resonant to the first red side-band of center-of-mass mode, so that the effects 
from the spectator motional modes can be neglected because of the very large off-resonant 
reason[15]. The Hamiltonian for this system is 

H = Ho + Hint, (1) 

n 

Ho ^u{a+a+ 1/2) +LoY,(^iz/^, (2) 

i=l 

n „ 

Hint = E 2 {^+^e*''^"^"^^e-^('^-'^)* + H.c] (3) 
1=1 

where u and a"*" (a) are the frequency and ladder operators of the center-of-mass mode. 
u is the energy difference between the ground state |0) and the long-lived metastable excited 
state 1 1) of each ion. For simplicity we further assume that the Rabi frequency and Lamb- 
Dicke parameter for each ion are same. If all ions are cooled under Lamb-Dicke limit, then we 
can expand the Hamiltonian up to the first order in 77. Transforming the Hamiltonian (3) to 
the interaction picture with respect to Hq and making use of rotating wave approximation, 
we have. 

Hint = -i^via^J- - aJ+) (4) 

where J± = Jx^iJy with = ^"^^ (Jia (a = x, y, z) being the three components of total 
spin operator of n ions. 

E.q.(4) is namely the so-called JCM Hamiltonian for n ions interacting with a laser beam 
(or beams) in an ion trap. Compared with the traditional JCM Hamiltonian, the total spin 
ladder operators of n ions replaced the spin ladder operators of one ion. When any one of 
the n ions is excited, the quantum number of Jz component of total spin operator increases 
by one, and in the meantime, the quantum number of center-of-mass mode decreases by 
one. Prom this mode, we would naturally like to look forward to some similar properties as 
shown in standard JCM. 
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In general case, searching for the solution of E.q.(4) is an open question. Below we will 
only discuss the case that including only two trapped ions. In this case the last term in 
E.q.(2) can be rewritten as cuJ^, and the quantum number corresponding to component 
only have three values: 1, 0, -l(the value "1" means both the two ions are excited, "0" means 
one of the two ions is excited, and means both the two ions are in ground state); Noting 
that in this system, the total excitations (including electron excitations and excitations of 
center-of-mass mode) is a conservation, thus we can arrange the eigenstates of Hamiltonian 
Hq in the order of total excitations, 

|1, n-2)", |0, n-1)", I -1, n)"; | 1, n - | 0, n)"+\ | -1, n + 1)"+^; • • • 

(5) 



Where the first numbers in each Dirac symbols represent the quantum numbers of compo- 
nent and the second ones represent the excitations of center-of-mass mode. The superscripts 
above each Dirac symbols represent the total excitations. 

In the Hilbert space constituted by above basis, the interaction Hamiltonian (4) is a 
diagonal matrix made up of 3 x 3 sub-matrices. Every sub-matrix represents a subspace 
corresponding to a definite excitation. In the subspace including n excitations, we can get 
the exact propagator for Hamiltonian (4) after a tedious deduction, 

^[(n-l)cos/3 + n] -^/^sin/3 Vgl) [_ ^ + i] 



U 



2n- 



n- 1 
2n- 



j sin P 



cos P 



^^2P[-cos/3+l] y^sin/? ^[ncos/3+(n-l)] 



(6) 



where (3 = ^/2{2n - l)rjVLt. 



III. COLLAPSES AND REVIVALS 

In order to study the phenomena of collapses and revivals of the two trapped ions, 
imitating traditional process, we assume that initially the internal states of the two ions 
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are both in the ground states and the external motional state of center-of-mass mode in a 
coherent state | a) . Thus the whole initial state of the system is, 

oo 

p(0) = q{m')q*{m") \ -1 m'){-l m" \ (7) 

m',m"=0 

with q{m) = exp{—^ \ a P);^] are the coefficients for a coherent state expanding in number 
state basis. 

The density operator of the internal states can be obtained by tracing over the vibrational 
motion of ccntcr-of-mass mode. With above equations, we can get the diagonal elements of 
the density-matrix of internal states for this two-ion system, 

(^) = f2^i^2i-^^^(^+^)Mn) (8) 

n=l ^ ^ 
oo 
n=l 
oo ^ 

= p(0) + J2 .^^_;^p (^cos/? + n - l)Mn) (10) 

n=l ^ 



with 



I ^ |2n 

p(n) = e-l°l^^ (11) 
n! 

is the probability distribution of phonon number for the coherent state of center-of-mass 
mode. 

According to E.q.(8)-(11), we plotted the time evolution of diagonal elements of density- 
matrix of the internal states as Fig.l, which describes the time-evolution of the occupation- 
probability of electronic energy- levels. We choose the initial average phonon number | a p= 
8, and the typical parameters 77 = 0.1 and 0, — 27: x 500 kHz (These typical parameters are 
used throughout this paper, below we will no longer narrate them repeatedly). Prom these 
curves, we can sec clearly collapses and revivals for the electronic state-occupation of the 
two trapped ions. It is worthwhile to point out that, as in the traditional standard JCM, the 
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effect of collapses and revivals in here also increases with the glowing of the average phonon 
number. However the average phonon number may not be very large in here, because we 
assumed that the ions are cooled to be under the Lamb-Dick limit. In order to observe 
collapses and revivals in the case of large number of average phonons, we must look for the 
solution of the system far from Lamb-Dick limit. This is a problem awaiting to be solved. 

Another possibility to obtain collapses and revivals for this system is letting the initial 
state of the ions in a superposition and the vibratic motion in a thermal state or vacuum 
state. This seems to be a more reasonable assumption, because the vibratic motion of the 
ions seems more to incline to a thermal state in general. For simplicity, we will suppose the 
motional state being in a vacuum, because in present conditions we can almost cool the ions 
to their vibratic ground state. Thus the whole initial state of this system is, 

1^(0)) = (a |1) + be^^^ |0) + ce^^^ |-1)) ® |0)„,, (12) 

with a^-F&^-Fc^ = 1. 

Similar to the deduction for coherent initial state, we have the diagonal elements of the 
density matrix for this kind of initial condition, 

pn = ^[2 + cos(V3ei)]' (13) 
Poo = y sin^Vs^t) + b^ cos^^t) (14) 

= ^a^[l - cos(V3^f)]2 + b'^sin'^iCt) + (15) 
9 

with ^ = \/2rjQ.. Numerical calculation indicated that only reveals clear collapse and 

revival. However, we may expect more clear collapses and revivals when the number of 
superposition-ion increased. Because the degree of coherence between atoms will increase 
with the addition of ions. The picture for this case is presented in fig. 2. 

In addition, we studied, by way of parenthesis, the behaviors of internal state of the two 
ions for initial conditions that the two ions being ground and the vibratic motion being in 
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number, thermal and squeezing states. Some results displayed in fig. 3. The amplitude of 
Rabi oscillation for the occupation of internal state |tt) is almost between and Iwhen 
the motional state be in a number state. But when the motional state is in a thermal or a 
squeezing state, the occupation-probability of |tt) is much smaller. 



IV. COHERENCE AND SQUEEZING 

In quantum optics, coherence and squeezing are two important aspects to describe prop- 
erties of a radiative field. Coherence implies, in some sense, the interference between differ- 
ent photons. And squeezing is a phase-dependent reduction of the noise of the electric field 
strength (or, of one of the field quadratures) below the vacuum level[16]. Now we borrow 
these ideas to describe the properties of vibratic modes of trapped ions. In this case, they 
represent the coherence between phonons and the squeezing of position or momentum of 
COM mode respectively. 

A quantitative measure for coherence can be defined by a fraction of the number of 
coherent phonons, 

2 

(16) 

to the total phonon number, 

oo 

(n(i)) = J]np„„(i). (17) 

n=0 

As for squeezing, we define the two quadratures of COM motion as, 

x = ^(a + a+) (18) 
p = ig{a — a"*") 

where g may be regarded as the structure constant of a harmonic oscillator. A quantitative 
condition for squeezing is that the normally ordered variance is under zero: 

(: (AO)^ :> < (19) 
8 



(n(t)) 



coh 



'n(n+l) 



it) 



n=0 



where O means x or p in E.q.(18). The normally ordered variance for x and p can be written 
as, 

(: {^xf :> = i^n{t)) + 2/ I V(n + l)(n + 2)i?e[pn(n+2)(t)] [ (20) 

f oo 



,n=0 



(: (Ap)^ :> = 2/ {n{t)) - 2/ i ^ + + 2)i?e[Pn(n+2) W] 

ln=0 



(21) 



.n=0 



Prom above two equations, we can see that in order to obtain squeezing, the phonon 
density-matrix must exhibit nonvanishing off-diagonal elements Pn(n+i) and /or Pn(n+2)- 
There may exist some different possibilities to generate them, e.g. initially (1) internal 
states be in a mono-state and vibratic motion be in a coherent state (or other states with 
coherence, say, binomial state [17].); (2) motional state be in a thermal or vacuum (means 
no coherence), but internal state be in coherent superposition; (3) internal state be in a 
superposition and, in the same time, vibratic motion be in a coherent state. Below we will 
still analysis the two different types of initial condition discussed in section 111. 

For the first case, e.g. initially the internal states of the two ions are both in ground and 
motional state of COM mode is in a coherent state. After a comphcate deduction, we get 
the density operator of vibratic motion, 

[(Z + l)(Z + 2) (m + l)(m + 2)]'/' 



Plmit) 



(2Z + 3)(2m + 3) 



1 - cos (\/2l + 3^t) 

(/ + l)(m + l 



1/2 



1 - cos(y2^et)J . (/ + 2) (m + 2) + ^^^^ ^ ^^^^^ ^ 
sin [^21 + l^t^ sin [y2m + l^i) q{l + l)g*(m + 1) + [{21 - 1) (2m - 1)]"^ x 
/ cos {\/2l - lit^ +1-1 [mcos (V2m - l^t) + m - l] q{l)q*{m) (22) 



where q{m) is shown in E.q.(7). 



It is easily imaginable that the average phonon number < n{t) > doesn't depend on (p, 
but increases with the augmentation of initial average number or < n(0) >. Given a 
definite initial phonon number r^, < n{t) > reveals also obvious collapses and revivals as 
time passes (see Fig.4). Moreover, this effects also become clear with the increasing of the 
initial phonon number. In fact, the collapses and revivals of phonon number relate to the 
collapses and revivals of the occupation of internal states. During the time that the number 
of phonons collapses, its number almost doesn't change. So the occupation of internal states 
also doesn't change, and collapses occur at this time. In order to compare with the collapses 
and revivals of the occupation of internal states in section III, we choose the same initial 
condition and parameters in Fig.4 as in Fig.l. 

Coherence 7 = |<a>|^/ < n{t) > is also independent on tp, but rely on initial phonon 
number < n(0) > and time t. Coherence 7 oscillates damply on time and reduces rapidly 
with the increasing of initial phonon number(see Fig.5). Thus if we want to keep the co- 
herence of vibratic motion for some long times, the initial average phonons must be very 
small. 

The squeezing of position and momentum depends on < n(0) >, (p and time t simultane- 
ously.. But from an inspection of (20)~(22), we can first easily make sure that the minimums 
of <: (a x)"^ :> and <: (A pY :> with respect to ip occur at <^ = 0. Then by use of nu- 
merical calculation, we find that <: (A p)^ :> increases rapidly with the growing of initial 
average number (see Fig. 6). When < n{0) >> 2.0, the momentum-quadrature doesn't exist 
squeezing. The minimum of <: (A p)^ :> with respect to initial average number occurs at 
about < n(0) >= 0.51. Given (p and < n{0) >, <: (A p)^ :> is an oscillating function on 
time. In order to see squeezing clearly, we plotted the time-evolution of <: (A pY :> in 
Fig. 7. As a comparison, we also plotted the time-evolution of <: (A x)"^ :> in the same 
initial condition. The maximum squeezing of momentum-quadrature in this type of initial 
condition is about <: (A p)^ :>= — 0.424g'^ which occurs at about (p — 0, < n(0) >= 0.51 
and t = 343/iS. Similarly, we can determine the minimum of <: (A x)^ :> and its location. 

Next, we consider the case that the initial state of the system is Eq.(12). In this case we 
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can deduce the density operator of motional state only having the following nonvanishing 
elements, 



1 r / ^ \ 1 2 
Pooit) = -a^ ^2 + cos (^VS^tj J + 6^ cos^ {^t) + c' 

Pii(^) = ^a^sin^ (y^Ct) + b'^ sin^ iCt) 



(23) 



(24) 



P22{t) = -a^ 



1 — cos 



poi(t) = pl^{t) = ;^«^sin (Vsit^ cos {it) e"^' + be sin {^t) e'^'^'-'^'^ 



(25) 



(26) 



Po2{t) = p*2o{t) = ^ac 



1 - COS (Vs^t 



(27) 



Pi2{t) = P*2i{t) = ^ab 



1-cos (VSit)] sin {it) e''^' 



(28) 



Similar to the above discuss, we can find average phonon number < n{t) > doesn't 
depend on ipi and ip2, and increases with the add of a, b, and also occurs collapse and 
revival on time. 

The coherence 7 reaches its maximum when cos(2(/7i — (P2) — 1- And this maximum 
increases with the enhance of c. In the limits that a = 0, 6 — > ( but b ^ ) and c — > 1, 
coherence 7 — > 1. This implies that in the case of very small phonons( "weak field" ), the 
vibratic motion approaches to a coherent state. In fact, in the above limits, the density- 
matrix elements (23)~(28) can be written as. 



Poo{t) = 1 - {n{t)) 



(29) 



Pii{t) = {n{t)) 



(30) 



Poi(^) = P*io{t) = e 



{n{t)) 



1/2 



(31) 
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and all other elements vanished completely ( Noting {n{t)) — 6^sin^ {^t) ). These are just 
the density-matrix elements of coherent state up to the terms of the order \a\ and |Q;|^(see 
the second reference of [8]). 

An inspection of (20)~(21) and (23)~(28), in combination with numerical calcula- 
tion suggests that for any (pi, when (p2 — 0, b — and c — 0.96, the squeezing of 
momentum-quadrature reaches its maximum.. This maximum squeezing oscillate with a 
nearly equi- amplitude on time. When t satisfying cos(-\/3C0 — the maximum squeez- 
ing <: (Ap)^ :>= —0.438^'^ for momentum-quadrature has been obtained (see Fig.8). In 
fact, in the condition that 6 = 0,a— >0,c— >1, and cos(\/3^t) = —1, the nonvanishing 
density-matrix elements in (23)~(28) are, 

poo{t) = 1 - ^ < n{t) > (32) 
P22{t) = ^< n{t) > (33) 

P02(t) = Ploit) = ;^e^^V< n{t)> (34) 

where < n{t) >= 16a^/9 with t determined by cos{V3^t) = —1. These are just the density- 
matrix elements of squeezing vacuum state (see the second reference of [8]). Thus we also 
obtained a "weak-field" squeezing state in another group of initial parameters. In the case 
that cos(\/3^t) 7^ —1, the momentum quadrature also has squeezing but not reaches its 
maximum (see Fig 7). 

It is worthwhile to point out that the coherence and squeezing discussed above includes 
the rapid time dependence, which, of course, is meaningless when slowly varying quanti- 
ties are of interest. But we can slowen the frequency of time-dependence by decreasing 
Rabi frequency Jl. Moreover, in current experimental conditions, dealing with problems in 
microsecond is completely feasible. 
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V. SUMMARY AND CONCLUSION 



We have studied the dynamical behaviors of two trapped ions interacting with lasers 
resonant to the first red side-band of center-of-mass mode. In small Lamb-Dick parameters 
and under rotating-wave approximation, the dynamics of the system can be described by 
JCM-like mode. An exact analytic solution for this type of JCM was presented. The results 
indicated obvious collapses and revivals for the occupation of internal states in two different 
types of initial condition. These collapses and revivals become more visible when the average 
phonon number increases. 

Coherence and squeezing for vibratic motion of center-of-mass mode were also discussed 
in these two different types of initial condition. For one type of initial condition that the 
internal states being grounds and vibratic motion being in a coherent state, coherence can 
maintain long times when the initial average phonon number is very small. A maximum 
squeezing 42.4% for momentum-quadrature was obtained in appropriate parameters. For 
another type of initial condition that the internal state being in a superposition and the 
vibratic motion being vacuum, apparent expressions for coherent state and squeezing state 
were acquired in suitable conditions. The maximum squeezing for momentum-quadrature 
in this case is 43.8%. 

We expect the results in this paper can bring some useful helps to other problems, espe- 
cially the quantum computation or quantum information processing. Because the research 
for the dynamical properties of two trapped ions locate completely in the field of trapped- ion 
quantum computation. We also expect this paper can give rise to some other interesting 
discuss for system of multiple ions, such as squeezing and coherence of vibrational state, 
collapses and revivals in large Lamb-Dick parameters or for systems containing more than 
two ions. 
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Fig.l. The occupations Pii(t), Poo{t), P-i-i(^) of the internal state are given as functions 
of time. Initially the motional state of COM mode be in a coherent state with average phonon 
number | a p= 8 and the two ions be in ground state. 

Fig. 2. The time evolution of p_i_i(t) for initially the internal state be in a superposition 
of E.q.(12) with a — b — c — and motional state be in a vacuum. 

Fig. 3. The time evolution of pii{t) for initially the internal states of the two ions are both 
in grounds and motional states are in number, thermal and squeezing state. In three cases 
the average phonon number (n(0)) = 3. (a) number state, (b) thermal state, (c) squeezing 
state. 

Fig. 4. The average phonon number < n{t) > is shown as a function of time for the same 
initial conditions and parameters as in Fig.l. 

Fig. 5. The coherence 7 is shown as a function of time and initial coherent phonon number 
for the same initial conditions and parameters as in Fig.l. 

Fig. 6. The evolution of <: (A p)^ :> as a function of the initial average number and time 
in given (p — 0. 

Fig. 7. <: (A xy :> (dot line), <: (A p)^ :> (full line) are shown as a function of time 
for initial parameters < n(0) >— 0.51, ip — 0. When t — 343/is, momentum quadrature has 
a maximum squeezing <: (A p)^ :>= —O.A^Ag^. 

Fig. 8. (a). <: (A xY :> (dot line), <: (A pf :> (full line) are shown as a function of 
time for initial parameters c = 0.96, h — 0. When t satisfies cos{-\/2>^t) — —1, momentum 
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quadrature has a maximum squeezing <: (A pY :>— — 0.438g'^. 
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